NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



•Mils available under NASA sponsoisWp 

In the int »>Tv in -ar’ wi-’s dis- 
semination of 

Progf3in intin ^nd wiUwut iiiDliiiy 
for any usa made thereot.' 


QUARTERLY PROGRESS REPORT I 


e t - 1 QtO.4 6 

tR-/4>3 75-p 


ELECTROMAGNETIC DEEP-PROBING (lOO-XCOO KMS) 

OF THE EARTH'S INTERIOR FROM ARTIFICIAL SATELLITES: 
CONSTRAINTS ON THE REGIONAL EMPLACEMENT OF CRUSTAL RESOURCES 

NAS 5-26138 


(E81-10046) ELECTROMAGNETIC DEEP-PROBING N81-1340b 

(100-1000 KMS) OF THE EARTH'S INTERIOR FROM 
ARTIFICIAL SATELLITES: CONSTRAINTS ON THE 

itEGIONAL EMPLACEMENT OF CRUSTAL RESOURCES UociaS 

quarterly Progress Report, 1 Jul. - 30 Sep. 03/43 00046 


John F. Hermance 

Department of Geological Sciences 
Brown University 
Providence, RI 02912 


Report Due Date: September 30, 1980 

Date of Submission: October 20, 1980 

Period Reported July 1, 1980-September 30, 1980 


1 


.^^eceived 

OCT i} , 1900 

SIS/ 902.6 
Type jc 

M- 0 0 7 


Statement of Work 


O bjective 

The objecLlve of this investigation is to evaluate the applicability of 
elecU'omiignel lo deep-aouniling experiments using natural sources in the 
magnetosplicre bv ineorporat ing Mngsat data with other geophysical data. 

Axi>roiU'|\ 

The Investigator shall piu'sue the abi'»ve objective through an analyala of 
Magsat satellite data, grovuid-based magnetic observations, appropriate 
ret'erevu'o flold modols, and other satellite data. 

The objective will be pursued by seeking the optimal combination of 
observations which lend i irst to a global, and then lo a regional, 
eharacterixat ion of the conductivity oi the Karth's upper mantle. 

Tasks 

The rollowlng tasks shall be performed hv the investigator in rulfllhnent 
of the above objective: 

a. Use data trom Magsat satellite to eonstrnin a long“period 
global "rospouso t unction" for the average Karth at low latitudes over 
a period ranging from n hours to 27 days, 

b. Sytvcltronlste tbo Magsat data with low-latitude ground-based 
observatory data to determine the vertical gradient of the respective 
magnetic field components. Use the vertical gradient of tl\o appropriate 
components to liulopendent Iv ascertain the separation of external and 
tutevnal field ivnitrihutions. 

c. Segregate the Magsat electromagnetic "response functions" 
according tt- live tectonic regime at tlu* I'arth's surlacc and evaluate 
svstomat Ic differences bet\'?een regions having lateral scale sixes on 
iho order of 1000 km or greater. 

d> Theoreiit al ly evaluate problems of resolution and Interpretation 
Involving electromagnetic induction by temporally and rpatial ly-varying 
magnotospltor Ic sources in a rv'tating Inhomogeucvnis Karth as observed at 
arblcrarv points in space. Use tlu'se theoretical studies to constrain 
the interpretation of Mngsat data as well as to propose further applica- 
tions of satol I llc-bascd clcciromagnevic decp-somuling experiments. 

e, Integrate the regional rosponse functions with other geophysical 
data in order to concirnln the joint tuierprctat Ion of comprcbcnslve 
physical models. 

f. Prepare and submit to NASA periodic progress reports and a 
detailed final report documenting the result. s of this investigation. 


FOCUS OF ACTIVITY 


PersotMtgl 

This? quai't:«‘r we have woblUsed personnel for this project ami hired 
a Ceophyalcal Onta Aneilyat (Michael Rossen) to assist in computer program 
development and the analysis of data. Mr. Rossen, a Physics graduate, 
will he supported in pan from tl<e NASA/NAOSAT program. 

Analysis of Cround Based Observatory Data (Ta sk B) 

Activity has focussed on developing a strategra for interfacing ground- 
based ci'servatory data with the satellite data. We have analyaed the 
Intltuduial do( endence of magnetic storm disturbances for periods available 
from the World Data Center (Idbb and 1967) in an effort to ascertain the 
suitability ol assuming a harmonic dependence for the source field. 

ihroblem! The dynamic source field seems to have a much more compll™ 
eatod spatial and temporal dependenee than global induction workers have 
hitherto assumed. 

Attempted jlolut ion: We are attempting to synthesise a dynamical 

field model at high sampling rates (I min. to 2.5 min.) which cein be e.aipared 
with the satellite data in the satellite frame of reference. 

Thepret ica 1 Jllode 1 tUmu 1 a t lpi\ ,CTa sk Jt) 

We are in the process of evaluating the coupling of lateral hotero'- 
genoities vs’lthliv the eartli Co.g. oceans) to source fields of finite 
dimensions. The vector components of the disturbance field are to be 
caleulnted at ground-level and at the satellite. The results of this 
caleulation will provide a basis for deciding the optimal use of ground- 
basod observations to a) minlmi«e noise at the satellite altitude; b) 
the feasihlUty of evaluating a globally averaged response function for 
a laterally heterogeneous earth; c) methods for regionalizing the re.sponse 
funetion of the earth to analyze lithospheric tectonics. 

Pro blem ; The theoretical analysis in three-dimensions of global 
slsed heterogeneities and arbitrary source-fields, easily over-runs the 
capability of most computer systems, particularly to model small scale 
features (the ocean-land interface) in ti>e global framework. 


a : 


Solu ti on ; Wo have reduced the problem to a two-dimensional model 
which seems to be valid at the depth of penetration and radius of curve- 
tui*e relations which pertain. This has to be evaluated! We ate beginning 
this analysis with the E-polarlsation mode and will be looking at some 
very simple sviurce field and geometrical relationships to gain insight 
into the problem. 

The Principal Inve.stigator conferred with .loe Cain when visiting 
henvor on other business and was appraised of certain problems of low- 
amplitude noise on the MACSAT data. Because of the relatively low level 
of magnetic disturbance during the MAUSAT mission wc are postponing the 
actual analysis of any MAUSAT data until the noise problem is mlnimixed. 

Filtering programs, developed by our group, for tlie analysis of 
space/tlme series have been forwarded to hr. Cain and are available for 
vise by other MAlhSAT invest igntors. Program listings are atenehed. 
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SUDFOtniNE FILT (X,N,DEL, FRg, PASS, HIP* SS ,XX # W) 
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DATA IN THE OBJECT 


X IS TUB ORIUINAL DATA 

THE SUruOUrXKE CCIJVERTS X (I) TO FILTERED 
N IS THE TOTAL NUSEEP OF DATA POINTS 
DEL IS THE SAHPLINO INTERVAL IN UNIT TTftE 
FRO IS THF HALF PO«ER (AMF) POINT IN CYCLES PER 
FILTER IS UEKERATED PY CONVOLVING ORIGINAL DATA 
LCE PASS OUTPUT IS THE SMCCTHED DATA 
HIGH PASS OUTPUT IS THE ORIGINAL DATA MINUS THE 
IADS IS THE SIGNAL AETER A LGP-PASS FILTER 
IIIFASS IS THE SIGNAL AFTFP A HIGH PASS FILTER 
TAILS CN THE ORIGIKAI DATA ABF CRFATFD PY REITERATING 
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UNIT 

WITH 


FIL0QO10 
PTL00020 
PROGRAFIL0Q030 
riL00040 
FILOOOSO 
FIL00060 
FIL00070 
FIL00080 
F1L00090 
EIL00100 
FILOOnO 
FIRST AND LAFIL00120 


TIME 

A GAUSSIAN 


SMOOTHFD DATA 


rCINTS. THESE TAXIS ARE SUrPFESSFD IN THE FILTER OUTPUT. FIL00130 
DIMENSION X(T) ,111PAKS(1) ,PASS(1) ,XX(1) ,W(1) FIL00140 

SIGHAsSORT( .69315 / (2. 0* (3. 14 1 6**2. 0) * ( FRO^'tB. 0) ) ) FIL00150 

TF*Sv'FT {f>.0*2.0t (SIUt!A**2-0) ) FILOQIsO 

«=TH/DFL FIL00170 

FiLOOiaO 

TP=AM*DEL FIL00190 

PlI 00200 

KIsNM F1L00210 

K«1*NtM+1 FIL00220 

KMM=tUH^M FIL00230 

FIL00240 

M«1sM + M^1 FILQt)25Q 

5UM=0.0 FIL00260 

DC 3 I*1,H FIL00270 

xxa)«x(i) piLOoaeo 

CCNTTNUE FIL00290 

PC 4 I-NFl,KKfl FILOOJOO 

XX(I)sX(N) FIL00310 

CCNTINUF. P1L00J20 

DC 5 1=1, N FIL00330 

K=IfK FIL00340 

XX(K)=X(I) FIL00350 

CCN'’‘INUS PI LO 03 60 

PC 6 K=1,MH1 FXL00370 

Cr.K-1 FILU03H0 

TC=C*nFL FIL00590 

U (K) =1.0/EXP(((TP-TC)t*2)/ (2.0* (S1GMA**2) ) ) FIL004 00 

CONTINUE F 11,004 10 

SUK=0.0 FIL00420 

rC 61 K=1,MM1 FIL00430 

SUN=W(K) +SilM FIT, 00440 

CONTINUE FIL004S0 

SUMA=0,0 FTL00460 

rC t3 K=1,MM1 FIL00470 

K (K)-R' (K)/SUH F1L00480 

SUMA=W (K) *SUMA FIL00490 

CONTINUE FILC0500 

PC B T=1,N FIL00510 

SUMA=0.0 FIL00S20 

no 7 K=1,MM1 FII.00530 

j=T+K-1 FIL00540 

snKA=SUMA* XX (J) (K) FIL00550 


>UE: MIT 


FOt^lFAH A 


♦♦♦ rcown University Coiputec Center 


7 CONTINUE rU00‘i60 

PASS(I)*SUNA moo«i70 

HIPATS -SUNA BILQO'^80 

8 CONTINUE FIL00590 

10 CONTINUE FII.UOuOO 

12 CONTINUE FliOOOlO 

11 TETUNN FtL00620 

end F1L00U20 


:itrz PAss*4 


FORTRAN A 


ni'own Univecssity Co«putf>c C«.*ntoc 


it 

c 

t 

c 


V. 

it 


c 

c 

r 


% 


V. 


!' 


. H 

\. 

V# 


V 




i,. 


V. 




SUBROtlTIlIK PASS4 (X, XHFh ,KRIUHT, KEFP, K, KMOD, DI, DTMOD , 
♦ PERIOD, SEL,LXTNn ,LErRS) 


ROUTINE BAND-PASS FILTERS BY CONVOLVING THE INPUT SERIFS WITH 
A DAMPED SINE WAVE. RPEOUEKCV RESPCNSE IS GAUSSIAN SHAPED, 
CEKTERFD AT THE GIVEN fPIOUENCY, WITH A NORMALIZED SELECTIVITY, 
GIVEN FY SEL . 

SUPfRCGRAMS RKCUIP’^O: 

NONE 


DIKERSIOK Xn» ,X»tW(1) ,KFIGHT(1) 


SUPRCUTINF PARAMETERS: 

X -ARRAY CCNTATKlhG l-CTI! INPUT AND OUTPUT SERIES OF EILTEK 

XNEW -TR.'UORAI'Y STOnAGF ARRAY (CONTAINS EXTENDED INPUT 
SERIES, SnOULC FE LAFUE ENOUGH TO ALLOW FOR THIS) 
WEIGHT-TEMPORARY STOFAGt ARRAY (CONTAINS FILTER FUNCTION, 

DIMENSION SHcni C IF GRFATFR THAN ( J *PBRIOD) /S EL*DT) ♦ 1 
NEFF -THE EFFICIFNCY PAl.AMETFF: THE FILTER OUTPUTS POINTS 
AT INTERVALS NEFF*i'T (NEFF IS USUALLY 1) 

N -NUMBER OR POINTS IR INPUT SERIFS 

NMOn -NUMBER OF POINTS IN OUTPUT SEFilES (NMOD= (N-1) /NFFE^I) 

DT -THE SAMFLING IMFRVM FBTWFEN POINTS IK INPUT SERIES 
DTMOD -IHF SAMPLING INTERVAL BETWEEN POINTS IN OUTPUT SERIFS 
PFFICD-CENTKR PElMCn (CCNVfcRTFC IMMEnlATELY TO FRQ- )/PER10C) 
SFL -SELECTIVITY OF THE FILTER: IF FRO(F) IS THE FRFOUENCY 
AT WHICH THE FILTER RFSTONSE IS DOWN TO 1/E OF ITS 
MAXIMUM VALUE, THFN SEL= (FRC-FRQ (E) ) /FLQ 
LXTND -PARAMETER CCNTRCLLi Nil TREATMENT OF THE ENPS OF X: 

-1 ZERO MEAN IS XSSDMFD FOR X, ZEROS ARE ADDED TO 
lACH END CUT TO 1/?. THE FILTER LENG'^H 
-2 X IS AVFRAGPr OVER ONE PFRICD INTO THE SERIES 
AND THE MIAN USFC TO FXTEKD THE INPUT SERIFS 
LSPRS -PARAMETER CCNTnCLLINv? TREATMENT OF OUTPUT SERItS: 

-0 EACH POINT WITHIN 1 FERICD OF EACH END OF THF 
SERIES IS MULTIPLIED BY ZERO TO MINIMIZE 
TRANSIENTS INlFOEUCED BY THF FILTER 
= 1 MO KODIFICATICN IS P’^RFORMED 


FIs3. mi5‘ju 

FRC= 1.0/PERIOD 

-> 

-> TC IS THE FORWARD LENGTH OF THE FILTER IN TIME 

-> CONSTANTS OF THF RAMPING TERM 

-> 

TC-2.0 

-> 

-> S IS THE UNMCRMAITZFl) SELtClIVITY 

-> 

S = SEI>FR'C*2.0*Pr 


PAS00010 
PAS00020 
PAS00030 
PAS00040 
PAS00050 
PAS00060 
PAS00070 
PASQOOnO 
PASOOORO 
PAS00100 
PAS00110 
PAS00120 
PAS001J0 
PASQ0140 
PAS0Q150 
PAS00160 
PAS00170 
PAS001B0 
PAS00190 
PAS00200 
PAS00210 
PAS00220 
PAS00230 
PAS00240 
PAS002P0 
PAS002fi0 
PAS00270 
PAS00280 
PAS00290 
PASOOJOO 
PAS00310 
PAS00J20 
PAS00330 
PASO 0.3 40 
PAS00350 
PAS00360 
PAS00370 
PAS00380 
PAS00390 
PAS00400 
PAS004 10 
PAS00420 
PAS00430 
PAS00440 
PAS00450 
PAS0046Q 
PAS00470 
PAS0Q4B0 
PAS00490 
PAS00500 
PAS00510 
PAS00520 
PAS00530 
PAS00S40 
PASOOS.SO 
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FILE; PASS4 


FonaaAK a 


♦♦♦ Brown University computer Center ♦** 


c > 

C > TFILT IS THE FOPHAPD LENGTH CF THE FILTER IN REAL TIME 

TF1LT»TC*4.0/S 

C— ”-> 

c > OMEGA IS THE CENTER FREQUENCY IN RADIANS PER UNIT TIME 

C > 


c 

c 

c- 

c 

c 


f 

r 

r*, 


C 

U 

C 


CMEGA=2.0*PI*FRQ 

> 

> MfILT IS THE HOMBEP OF SAMPLING PITS IN 1/2 CF THE FILTER 

> TFILT IS ADJUSTED LENGTH CF THE SERIES 

> KFILI IS TOTAL NUHDFr. CF SAMPLING BITS IN THE FILTER 

HFILT-=TFILT/DT 

TFILT^MFILT 

IFILTsTFILT^DT 

NFILT=2*MFILT«-1 

AEFF=NEPF 

•--> 

> THE FILTER FUNCTICN IS GENERATED 

> 


DC 10 1=1,NFILT 
TaI-1 

T=T*DT-TFILT 

WEIGHT (I) =2.0*S* (SORT (PI) ) *CC5 (CKEGA^T) /EXP (S*5*T*T/4.0) 
10 CONTINUE 
> 

> THE ENDS ARE EXTENDED 

> 


IF(LXTND.EQ.2) GO TO 11 

XEND=0.0 

XEEGIN=0.0 

IF (LXTND.EQ. 1) GO TO 21 
11 KEX=1.0/FnQ/DT 
KEX=KBX+1 
SUE=0.0 
AFX=KEX 

DC 100 1=1, KEX 
SUM=X (I) ♦SUM 
100 CONTINUE 

XBFGIN=SUM/AEX 
KKEX=N-KEX* 1 
SUM = 0,0 


DC 110 I=KKEX,N 
SUM=X(I) ♦SUM 
no CCKIINUE 

XEND=SUM/AEX 

> 

> NEW IS THE NUMBER Of POINTS IN THE NEW SERIES 


21 NEW=N*2*MFILT 
NEND=MFILT+N 
DC 20 1=1, MFILT 
IDEGTN=I 
IENB=I+NENn 


p'V'iS^ O'.'ALrrV 


PAS00560 

PAS00570 

PAS00580 

PAS00090 

PAS00600 

FAS00610 

PAS00620 

PAS00630 

PAS00640 

PAS00650 

PAS00660 

PAS00670 

PAS006d0 

PAS00690 

PAS00700 

PAS00710 

PAS00720 

PAS00730 

PAS00740 

PAS00750 

PAS00760 

PAS00770 

PAS007fl0 

PAS00790 

PAS00800 

PASOOaiO 

PAS00820 

PAS00830 

PAS00340 

PAS00850 

PAS00860 

PAS00870 

PAS00B8O 

PAS00890 

PAS00900 

PAS00910 

PAS00920 

PAS00930 

PAS00940 

PAS00950 

PAS00960 

PAS00970 

PAS00980 

PAS00y90 

PAS01000 

PAS01010 

PAS01020 

PAS01030 

PAS01040 

PAS01050 

PAS01060 

PAS01070 

PAS01080 

PAS01OS0 

PAS01100 


r> f > r, n o o o o o o n 


PILE: PASS4 


FO&TBAM A 


♦♦♦ Brown Universifcy Conputtc center 


XKEB (IEIGIM)»XBIGIS 
XKEH (lEHD) *XEKD 
20 CCNTINOE 

— > 

> THE OBIGINAL SERIES IS ADDED TO THE fllDDLE OF THE NEW SERIES 

> KNEW IS THE EXTENDED SERIES 

> 

DC 30 1*1, N 
O^I*HPILT 
XNEW(J)«X(I) 

30 CONTINUE 

> THE EXTENDED SERIES IS CONVOLVED WITH THE 

> FILTER IMPULSE RESPCNSE FONCTICN 

KHOD=0 

I«CD=0 

DO 40 I=1,N»NEFF 
SUM=0.0 

DC *>0 JF1LT=1,NFIIT 
K*I4JFILI-1 

SU«==SUH4XKEW (K) *WEIGHT(JFILT) 

50 CONTINUE 
IMCD=IH0D4l 
NECD=iaCD 

X (IMOD) sDT^SUM/ (2.0*PI) 

40 CONTINUE 

IF(LSPRS.EQ. 1) GO TO 31 

> LTRANS IS THE LENGTH CF THE TRANSIENT 

LTRANS-PF.RIOD/ (DT^AEFF) 

LDEGIN=LTRANS4l 
IEND=HHOD-LTRANS 
LEKDl-LEKDfl 
DO 60 1=1, LTRANS 
X (I) =0.00001 
60 CONTINUE 

DC 80 I=1,LEKD1,N«0D 
X (I) =0.00001 
80 CONTINUE 
n CTKOD=Dr*NEFF 
LEIURM 
END 


PAS01110 
PAS01120 
PAS0H30 
PAS01 140 
PAS01150 
PASO 1 160 
PAS01170 
PAS01180 
PAS01190 
PAS01200 
PAS01210 
PAS01220 
PAS01230 
PAS01240 
PAS01250 
PAS01260 
PAS01270 
PAS01280 
PAS01290 
PAS01300 
PAS01310 
PAS01320 
PAS01330 
PAS01340 
PAS01350 
PAS01360 
PAS01370 
PAS01300 
PAS01390 
PAS01400 
PAS01410 
PAS01420 
PAS01430 
PAS01440 
PAS0 1450 
PAS01460 
PAS01470 
PAS014C0 
PAS01490 
PAS01500 
PAS01510 
PAS01520 
PAS01530 
PAS01540 


ri r> «■> f V r» f> r> <■> o n n» o o f> rj o o r> n o o r» r» o <tK n 'i o <s- O' O' o r> o r?' r» p 


I IMIltTS 


SUB 8 00 TINE WIM® 




I IS THE ARRAY eORTAZNTiNMS lEllR 

ALSO THE FtlTRERD SSBXE§ t§ 
t If A IBB-PORAIiV S-TagAtl ARRAY, ftS 'St f |I1 

THE 8AIN PR'OBRAS) Mf-tT AXL®| 'fiS'l 

BY A. RERIOB LE«'B 'OR EA'CS ®iS fO 'liSii 

sti€l fKRii&i) 
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LX IS lES SUillB Of PGIR^a IN !• 

lY IS IBS NiMSIS ©? POIHtS CONfAINSB 18 X 
itif IS tBi SAigLINS IRfSBYAl, 

PIRIOB IS Till CENTER PBRIOS Of fHl BAB®-*FASS flltm®, 

S'tL IS Tit SiLiCTIVIft OF TiS fllftS, It IS ftfS'SBl 

AS (F2-K1) /FO MESS fO IS Tit 0I8TS1 OF TBB BA81{ 

FI AND F2 ARS THE FHSQUHNCItS T© lACi SIDE OP F© 

WiStE Tit AiPLITiDB RBSfONSE OP TIM tllTBt JS ItiAL 
TO 1/2. TiS RESPONSS AT f© IS 1. 

KE8D: sO, THE X ARRAY IS FILTERED ASSIBINS THAT Tit VALUSS 
ABE ZiSO ODTSX0H TSS ARRAY. TBS ENDS ABB NOT 
SXTBNDED. 

A CCSINS TAPER IS APPLIED T© Til INDS OP TiS INPif 
SlSIiS. THE ENDS ARB NOT EXTINDID* 

^2, THE ENDS OP THE INPUT SEBIBS ABE EXT8ND1.D BY OHB 
CENTER PtaiOD AND A COSISS TAPBE IS APPLIED TO TIS 
EXTENSION. 

OilllNTs THE FILTHS DOBS NOT HAVE 2SRC PiASE SlIFT, SO Til DATA 
IS FILTERED POSWARD IN TINS AND TiEN THi FILTSR IS 
APPLIJD IN Til RBVISSI DiRlOtlON TO OPPSET THE PiASl SiltT. 

SUBPlOOitAMS ESQUIRED: 

NGHB 


DIMENSION X(LIC),y(LY) 

CO MP tl X SO , SOC tZi, GlIE , C 
INTBOSR EXTEND, TAPEl 
PO-1,0/PBEIOD 

S ELICT = SEItSQRT ( 2. ) ♦DELT*FC* 1 . 13 94 
C=(G.0,1-Q) 

ONF,= (1.0,0.0) 

Pl=3. 1415926 
F.X1END=P 

20= ( 1 . 0 ♦ SSL ECT) ♦CIX® (-2. 0*P1TC* F0*D1LT) 

2CC=CCKJG(iO) 

Z1=CEXP (-2. 0*PI*C*PG*DBLT) 

iO= (CABS (21-ZO) ♦CABS (tl-tOG) ) / (CBiS {Z1--GNI) *CABS (Zl^ONl) ) 
AG=HO/BEAL(ZC^aOC) 

A2 = A0 




PHI00070 

P0JOOD80 

P1I00090 

PHI00100 

PIJIQO110 

PHIQ012O 

Pni0Ql3O 

PHI0O14O 

PlISOlSO 

PiX0O16O 

paiooiTo 

PHIQD180 
PU100190 
PHI00200 
Pta0D21O 
PHI00220 
PIIQ0235 
POX00240 
PHIOO250 
PHI00260 
PHI00270 
PHI002RO 
PH10029© 
PHI 00 3 00 

PH100310 

PHI60320 

PHT00390 

PHI00400 

PHI00410 

PUI00420 

PHIO'04.30 

PH100440 

PH100450 

PHI00460 

PHI00470 

PHI00480 

PH100490 

PHI00500 


tllEj PHllTR FORTH A 


»*» ncoi#ii University Computer Center **• 


P1»(HEAL(Z0f70C) )/ (FEAi (7C*ZCC) ) 

P2*1.0/(REAI. (Z0*20C)) 

TAPFR * 1 
IXKEWLX 

IF(KEND-EQ-O) GO TO 35 
IF(KEKD.EC. 1) GC TO 25 
EXTBKD*i?ERIOD/UELT*>1 

^ — > EXTEND IHH ENDS OF Tilt INPUT TIfJE SEflIfS 

DC 10 ID*1,LX 

X(r.X*FXTEKD-IL^1)*X(tX-IL>1) 

10 CGNTINIIE 

DC 20 lEal, EXTEND 
X(1E)»X(EXTEND^1) 

X (tX^EXTF.KD4lE)*X(LXfEXTFKD) 

'0 CONTINUE 

lXKEK-LXf2*EXTENC 
.5 CONTINUE 

; > APPLY A COGINE TAPED TO EITHER THE INPUT OK EXTENDED 

> Tinr SERIES 

TAPER* PEraCT/(3.0*DtlT) 41 
ro 30 IT=1, TAPER 

THETA= (PI/2.0) ♦ (TAI-r:P-IT4l) /TAPER 
X (IT)=X(TAPFI'4l)*COa (IHFTA) 

X(LXNEW-IT4l)*X (LXNEX-TAPRR) *C0S(TI1RTA) 

'0 CCKTINHE 

■5 CONTINUE 

• -> FILTER IN THE DIRECTION OF INCFEASING TIME 

X (1)=-AO*X (1) 

Y(2)=-AC*X(2) 4EUY (1) 

DC HO IL=3,LXNEW 

Y (IL)=-AC*X(IL) 4A2*X (lL-2) 4P14>Y (IL-1) -B2^Y (IL-2) 

10 CONTINUE 

DC 50 IT-1, TArER 

T!1ETA= (PI/2.0) ♦ (TAPf.P-IT4l) /TAPER 

Y (TT) =Y (TAPERfI) ♦COS (THETA) 

y; (LXKEW-IT4l)=Y (LXNEW-TAPEi:) ♦CCS (THETA) 

.0 CONTINUE 

> 

• > ULTER IN THE DIKtCIiCK CF DECREASING TIMS 

> 

X (LXNEW) =-AO^Y (LXNEK) 

X (1XNEW-1)=-A0^Y (LXNEH-1) 4D1 *X (I.XNEH) 

DC 60 IL*3,LXNBfc 

X (LXNEW-IL4l)=-A0^Y (LXNEH-IL+1) 4A2^Y (LXNEH-II.43) 

♦ ♦D1»X (LXNEH-IL42) -B2^X (LXNEW-IL43) 

0 CCKTINUE 

DC 70 IL*1,LX 
X (IL) =X (EXTFND+IL) 

’0 CCKIINUE 


PHI00510 

PHI00520 

PHI00530 
PHI00540 
PHI00550 
PH I 00 560 
PHI00570 
PH100500 
PH I 005 90 
PHI00600 
PHI00610 
PHI00620 
PHI00630 
P11I00640 
PHI00650 
PHI00660 
PUI00670 
PHI00680 
PHI00690 
PHI00700 
PHI00710 
PHI00720 
PHI0073Q 
PHI00740 
PHI00750 
PHI00760 
PHI00770 
PHI00780 
rHl00790 
PHI00800 
PHI00810 
PHl00fJ20 
PIII00830 
PHI00840 
PHI00850 
PHI 00 860 
PHI00870 
pHiooeeo 

PH1C0890 

PHI00900 

PHI00910 

PHI00920 

PHI00930 

PHI00940 

PHI00950 

PHI00960 

PHI00970 

PIII00980 

PHI00990 

PHioiono 

PlllOlOlO 

PHI01020 

PHI01030 

PHI01040 


fILl: 

« 


PHIIIB POBIBAU A 


*♦* Bcown University Cosputer Center *** 


BEIUBH 

END 


PHI 01 050 
PHIQ1060 


r, n n r. r'J rt n r. 



] 

i 


yiLf: SMOwlH FOKTHAK A Crown University Computer Center ♦♦♦ | 


c snoooolo 

SOFBOUTIKE SHOOTII(X,K,«) S«000020 | 

C SHO00U30 

sMuooaao I 

PPOGBAM 10 SHCOTH DAIA USIKG EQUALLY WEIGHTED RUKN1K6 AVERAGF SEOOOOSO I 

SKOOOi)60 

X IS THE IKPUT DATA SECaES. SH000070 

H IS THE NUHDEk CE INFUI DATA POINTS. SMO00J80 

M IS IIIK HIEEEA OF PCIKTS IN THE AVERAGING INTERVAL. SM000090 

snoocioo ■ 

SHC00110 

niMENSlCN X(1) ,SflUTH(3000) Sf5000120 

A*« SHO00130 

1EEG1N»(}1*1) /2 SHC00140 

IEKn*N-IPEGINf 1 SMO00150 

L««(?!-l)/2 SMC00160 

DC 10 I=IfiEGlK,lEND SMC00170 

SUMsO.O SHCOOiaO 

DC 20 LSMUTH«1,tt SM000190 

I,*LSf.UTH-l SMO00200 

3»I*L-LM SHO00210 

SUK=SUH+X(J) SH000220 

20 CONTINUE SM000230 

SBUTH (I)=SUK/A SB000240 

10 CONTINUE SMC00250 

r SHO00J60 

; THIS STEP OVEREMPHASIZES SMOOTHING AT ENDPOINTS. SMCOOITO 

DC 30 1^1, N SMO00280 

IF {I.LT.IBEGIN) SMUTH (I)-SMUTU (IBEGIN) SM000290 

IF (I.GT. lEND) SMUTII (I) =SKUTII ( IFND) SMOUOJOO 

30 CONTINUE SMO00J10 

C SM000320 

C THIS STEP DESTRGXS THE ORIGINAL DATA SERIES, REPLACING IT SMO00330 

C WITH THE SMOOTHED SERIES. SM000140 

DC 40 1=1, N SKCOOiOO 

X (I) =SKfITH(I) SM000360 

40 CCNTIHOE SMC00370 

RETURN SM000380 

END SMO00390 

SUBROUTINE SPLIT (N , ISIT) SHO00400 


t 



